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We study the damping of the collective excitations of Bose-Einstein condensates in a harmonic 
trap potential loaded in an optical lattice. In the presence of a confining potential the system 
is non-homogeneous and the collective excitations are characterized by a set of discrete confined 
phonon-like excitations. We derive a general convenient analytical description for the damping 
rate, which takes into account, the trapping potential and the optical lattice, for the Landau and 
Beliaev processes at any temperature, T. At high temperature or weak spatial confinement, we 
show that both mechanisms display linear dependence on T. In the quantum limit, we found 
that the Landau damping is exponentially suppressed at low temperatures and the total damping 
is independent of T. Our theoretical predictions for the damping rate under thermal regime is 
in completely correspondence with the experimental values reported for ID condensate of sodium 
atoms. We show that the laser intensity can tune the collision process, allowing a resonant effect for 
the condensate lifetime. Also, we study the influence of the attractive or repulsive non-linear terms 
on the decay rate of the collective excitations. A general expression of the renormalized Goldstone 
frequency has been obtained as a function of the ID non-linear self-interaction parameter, laser 
intensity and temperature. 

PACS numbers: 03.75.Be, 03.75.Lm, 05.45.Yv 


I. INTRODUCTION 


The damping process plays a crucial role in the dy¬ 
namic of Bose-Einstein condensates (BEC). Phenom¬ 
ena such as superfluid phase transitions™ 7 Josephson 
effect quantized vortices, 10 ^ Mott insulator transi¬ 
tion,^ among others,^ are limited by the finite lifetime 
of the collective excitations of the condensed atoms, i.e. 
by the damping mechanisms and their dependence with 
the temperature. After the experiment al co nfirmation 
that the collective excitations are damped J2 MU 5 U 6 1 the kg. 
havior of the decay rates have been of mayor interest in 
the physics of BEC. Thus, several studies and theoretical 
calculations of the damping of e xcitat ions have been per- 
for med in three-dimensional (3D)R-Ll1122J, two-dimensional 
(2D)fll22if2l as we n as one-dimensional(lD) systems!^ 27 ! 
For a better understanding of the damping process, it 
becomes necessary to consider the contribution of the 
parabolic confining potential. The behavior and charac¬ 
teristics of the decay rate and its dependence with tem¬ 
perature, differer radically if we are dealing or not with 
homogeneous systems. The main assumption for homo¬ 
geneous system is to consider the condensate density con¬ 
stant in the all space. 

In the present work we are dealing with a micro¬ 
scopic theory for the damping rate of collective oscilla¬ 
tions, specifically for a quasi-one-dimensional (ID) con¬ 
densate confined to a parabolic harmonic trap potential 
and loaded into an optical lattice. The existence of a non- 


negligible confining external potential breaks the invari¬ 
ance symmetry, which leads to the damping rate showing 
a different qualitative behavior in comparison with pre¬ 
vious formalism, where the condensate is tackled as a 
homogeneous system (see for example Refs. [23 and ESj 
and references there in). The influence of both exter¬ 
nal interactions -the trap potential and laser intensity- 
must provide a physically richer scenario for the decay of 
collective excitations. The knowledge of excited states or 
Goldstone modes enables the characterization of the con¬ 
densate dynamics in a general framework. In the case we 
are in presence of spatially non-homogeneous BEC sys¬ 
tem, the label spacing of the discrete spectrum (confined 
phonon-like modes) and also, the nature and symmetry 
of wavefunction of the collective modes, are required for 
the calculation of the collision scattering process.^ 
Cigar-shaped traps can be considered as quasi-one- 
dimensional systems. We select such platform of a con¬ 
densate loaded simultaneously into a ID harmonic poten¬ 
tial and an optical lattice to characterize the phenomenon 
of damping and tackle the problem analytically. Within 
the framework of mean field theory, the physical charac¬ 
teristics of a BEC in such trapping profile are ruled by 
the time dependent non-linear Gross-Pitaevskii equation 
(GPE) [MUD j n an external potential 
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where m is the alkaline atom mass, Vl the laser inten- 
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order, the collision term is described by the interaction 
between three interacting phonon modes, giving rise to 
the vertices in the self-energy diagrams shown in Fig. [ljo) . 
It becomes clear that the last diagram in Fig E>) does 
not contribute to the self-energy interaction at T = 0 
K, since a thermal excited mode u>i must be present in 
the system (see Eq. ([3]) below). We must recall that the 
thermal cloud in the present theory is assumed in thermal 
equilibrium. 

Figure So presents the leading diagrams contributing 
to the self-energy 7r p . Accordingly, the complex frequency 
correction can be factorized into two main processes as 
7 r p = ACjl + A Cjb and, in the Hartree-Fock-Bogoliubov 
approximation we obtain that 23 


2' 7r <?i 

H 2 ^ 


2 (fi ~ fj) \Aij\ 2 

ujp + uji — ujj + is 


(1 + fi + fj) \Bjj \ 2 
oj p — uji — uij + is 
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FIG. 1. a) Diagrammatic equation for casual Green function 
G p of phonon modes in the condensate, b) Vertices involved 
in the collision phonon-phonon cubic interaction, c) Feynman 
diagrams contributing to phonon self-energy operator n p for 
Beliaev and Landau processes. 


where gi is the coupling constant, A l3 (Bjj) represents 
the matrix elements for the Landau (Beliaev) process 
oj p +uJi —> u)j (ujp -A LUi+LOj). In Eq. (31 the sum Y2 takes 


into account all possible virtual transitions | wf) and | ujj ),, 
contributing to the decay rate, while the term 1 + fi + fj 
gives us the Bose-Einstein statistical factor of the phonon 
|w p ) decaying into two confined phonon modes (first dia¬ 
gram in Fig. |TJc)) |uij) and | uij) , and /, — /)■ corresponds to 
the thermal correction for the annihilation and creation 
of phonons with frequencies uii and ujj (second diagram 
in Fig.[iJ)), respectively. 


sity, d its laser wavelength, and wo the frequency of the 
harmonic trap. 

II. THEORETICAL BACKGROUND 

In the framework of the Green function formalism, the 
spectrum of the excited states is obtained by the poles 
of the dressed Green function G p . The solution of the 
Dyson equation, shown diagrammatically in Fig. EO, is 
the renormalized Green function G p given by 

G p = G 0p — 7 T p . (2) 

In absence of interaction G p -A Gq p = [ui — oj p + zer] — 1 , 
s > 0, uj p is the eigenfrequency of the excited state, and 
7 t p is the self-energy contribution. The solution of Eq. (pi) 
leads to the complex frequency u> = oj p - t-7r p . Here, Re{jKp) 
represents the renormalized contribution to the eigenfre¬ 
quency uip, while Im (—7r p ) corresponds to the damping 
rate. 

We assume that the decay processes are associated to 
the collision between confined phonon states. In first 


A. Bogoliubov type excitations 

Bogoliubov type excitations can be searched by apply¬ 
ing a small deviation from the GPE stationary solutions 
I'i/’o) exp(— ifit/h), i.e. 

l'F(t)) = exp(—int/h) [\ipo) + |u) exp(— icot) 

+ |t>*) exp(iwi)] , (4) 

with g being the chemical potential. By lineariz¬ 
ing the time-dependent non-linear GPE, we obtain the 
Bogoliubov-de Gennes equations (B-dGE) for the eigen- 
frequencies oj and amplitudes |it) and |u). As stated 
above, if the harmonic trap potential is switched off, 
we are in presence of the homogeneous case, where the 
phonon wavevector q = q x e x is a good quantum num¬ 
ber and with excited frequency co(q x ) depending on the 
phonon wavevector. Thus, the Bogoliubov’s low-lying ex¬ 
citation spectrum shows a linear phonon dispersion in q x . 
On the other hand when the condensate is loaded into a 
harmonic potential, wo / 0 in Eq. ([I]), the system be¬ 
comes inhomogeneous and the wavevector q is no longer a 
good quantum number. In such case the B-dGE provides 
a set of discrete excited state frequencies uj p (p = 1, 2,...). 
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By considering the periodic potential in Q and the 
non-linear term, g\ \4>o\ 2 1 in the GPE as a perturbation 
with respect to the harmonic trap potential, ^mw^a; 2 , 
phonon frequencies w p up to the second order in gi are 
given bySS! 
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with A = g\N / (IqTujJq), N the number of atoms, Iq = 
^h/mujQ defining the characteristic unit length, a = 
2 tt1 0 /cI and Vo = Vl/Swo- In addition L t (z), r(z), Ei(z), 
Chi(z), and C are the Laguerre polynomials, the gamma 
function, the exponential integral, the cosine hyperbolic 
integral and the Euler’s constant, respectively. The pa¬ 
rameters 7 p , S p and p p are reported elsewhere.^ 



FIG. 2. (Color online) Dependence of the label spacing 
Ap = (w p — in — ujj ) /wo and A]*’ p) = (uj — uh — u> p ) /w 0 
(see Eq. §) on the reduced laser intensity Vo for Beliaev and 
Landau damping rates, respectively. Critical values, Vg P ’ 1 ’^, 
where Aapproaches zero, are shown by arrows. 


B. Symmetry of the excited states 

Owing to the inversion symmetry, the space of solu¬ 
tions can be decoupled into two independent subspaces, 
O and £ for p = 1,3,.... and p = 2,4,.... modes, re¬ 
spectively. Hence, the components \u p ) and |u p ) are ex¬ 
panded over the complete set ID oscillator wave func¬ 
tions {(j> 2 p+i} or {(f> 2 p } for the odd, O and even, £ 
Hilbert subspaces. The normalized eigenvectors |<h p )' = 
[|Q , |Up)] , up to first order in A and Vo, can be cast 

a#a 
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/ IJ. \ I A/p m Vo g p m) | ^ N \ 

\<Pp) + L -d?-1 

m^p 2 (p-m) 
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2 \(p-m)/2 

exp (-a 2 ) L lp ~ ml ( 2 a 2 ) , ( 8 ) 


L^(2a 2 ) the Associate Laguerre polynomials, h = (p 
m — \p — 7?i|)/2 and m + p is an even number. 


The parity of the function |$ p ) is linked to the in¬ 
dex p, if p is even or odd the eigenstate |$ p ) is sym¬ 
metric or antisymmetric. The decay process of a cer¬ 
tain phonon p is restricted by the symmetry property of 
the matrix elements in Eq. (13]). The amplitudes Aij(p) 
and Bij ( p ) impose a parity selection rule for the involved 
states |$p), | Ui) and |uj). As shown in Eqs 


( 111 , for a symmetric (antisymmetric) state 


• © 

\*p), 


and 
the 

amplitudes |uj) and |uj) must fulfill the parity condition 
i+j = even (odd) number, therefore limiting the possible 
number process for Beliaev, w p —> w, + ui,j. and Landau, 
Wp + Ui —> uij decay rates. Besides the symmetry of 
the matrix elements Aij(jp) and iy(p), for certain eigen¬ 
mode |<hp) with frequency w p , a key role in the damping 
process is ruled by the label spacing between the Bogoli- 


ubov collective oscillations A 


(m) 


= (Wp - Wi - OJj) /w 0 
and A/’ p) = (wj — Wj — w p ) /wo- Fixing the frequency 
Wp, all allowed combinations of w* and Uj approaching 


A p ’^ or A/ p) to zero value, lead to resonant transi¬ 
tions for the Beliaev or Landau damping processes. This 
effect is shown in Fig. [2] where the influence of the in¬ 
tensity Vo on some level spacing A p ’^ (panel a)) and 
A/’ p) (panel b)) are displayed. For calculations we fixed 
d/lo = 0.25 and A = 2. From the figure, it can be seen 
that the laser intensity can be used as a external parame¬ 
ter to tune particular damping process, i.e. we are able to 
reach certain critical values Vq P ’ 1 ’^ (see Fig.^a)), where 
A p "' ]> (V 0 ) = 0. This is a direct consequence of the fact 
that the Bogoliubov-type collective excitation energies, 
as a function of A and Vo, are not equidistant. Notice 
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that, for a certain state p , the Landau mechanism allows 
more combinations fulfilling the condition A^ l ’ p \Vo) = 0. 

III. DECAY RATE 

In the sequel, we consider that the damping is orig¬ 
inated by a collision process, and in first-order approx¬ 
imation, it is described by the interaction of the three 
confined phonons, giving r ise to a cub ic interaction in 
the bare phonon amplitude! 19 31 22 l 31 l 32 l This mechanism 
is represented by the vertex diagrams of the self-energy 
part shown in Fig. E>)- 


A. Beliaev damping 

In first-order loop approximation, the Beliaev mecha¬ 
nism arises from the collision of three particles where one 
phonon with frequency io p is annihilated decaying into 
two confined excitations uii and ujj. Therefore, the al¬ 
lowed processes for the confined modes oj p are those with 
p = 2,3,.... Following Feynman diagrams of Fig. fib) 
and using the eigenfunction amplitudes given in Eq . (|6|, 
we find that the decay amplitude B^ can be cast aip22IJ 
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FIG. 3. (Color online) Reduced Beliaev damping y^ f° r 
the Goldstone modes p=2,...,8 versus the dimensionless self¬ 
interaction parameter A at laser intensity Vo = 0 and T =0 
K. 

Bij = J dxipo [■Up (u*u* + u*v* + v*u*) + 

v p (u*v*+v*v*+v*u*)] . (9) 



FIG. 4. (Color online) Influence of the laser intensity Vo on 
the reduced Beliaev damping y^/yAb Resonant peaks are 
related to the zeros of the label spacing Ap ,jl (Vo) in Eq. ([ 3 J . 


Thus, for the Beliaev damping rate we obtain 

7^ ) = ^|A|H B) (A,y 0 ) > ( 10 ) 

with y(°) = 4:Trgi/(loh) and M.p B \A,Vo) being defined 
in the Appendix A . Notice that Beliaev mechanism is 
forbidden if e —» 0 in Eq. ([ 3 ] ). The energy conservation 
limits the real phonon transitions ui p —» w* + ujj. 

Figure [3] displays the behavior of the Beliaev damp¬ 
ing rate, 7 ^, in units of y*- 0 ^, as a function of the re¬ 
duced parameter A for the first seven allowed confined 
modes p = 2,...,8. In this calculation we used T = 0 
and laser intensity Vl = 0. For small values of A, all 
the normalized damping seen in Fig. [3] presents a lin¬ 
ear behavior, while for increasing values of |A|, the func¬ 
tion 7 ^/ 7 ^ behaves non-monotonically and reaching a 
maximum. For a given excited state |4> p ) , the position of 
the maximum is not symmetric with respect of the type 
of non-linear interaction (repulsive g± > 0 or attractive 
<71 < 0). In Fig. UJ we show the dependence of 7 ^ on 
the dimensionless laser intensity Vo, and calculated for 
A = 2, d/lo = 0.25 and T = 0 K. It can be observed that 
the excited states p = 4, 5 and 6 .show sharp peaks at cer¬ 
tain values of Vo These features are linked to the zeros 
of the frequency label spacing Ap"' J ' > (Vq) as represented 
in Fig. [2j while the number of transitions uj p —> ay + a jj 
and the strength of the matrix elements, B r:j , dictate the 
relative intensity of the peaks. 
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B. Landau damping 

Here, in the damping process a phonon mode with fre¬ 
quency and a thermal excitation cu, are annihilated 
and confined phonon is created. Thus, the Landau mech¬ 
anism is a thermal process at finite temperature. In 
present case the vertices phonon-phonon interaction (see 
Fig#)), conduces to the amplitude probability) 19 l 23 l 

Aij = J dxcj) o [• u p ( UiU * + ViV* + ViU*) + 

V p ( UiU* + ViV* + ViUj)] . ( 11 ) 

Thus, we have 

7 i p) =7 (0) |A|H i) (A,y 0 ) , (12) 

where A4p L) (A, Vo) is defined in the Appendix B. Hence, 
the total damping can be cast as 

7 (p) = 7 (0) |A| ( M^ + \m^ . (13) 

Figure [H] presents the total damping 7^ (solid lines) 



Beliaev damping, the first excited state p = 1 is forbidden 
at any temperature, while for T / OK this mode becomes 
allowed for the Landau process. 



FIG. 6 . (Color online) Influence of the laser intensity Vo on 
the total damping 7 ^ = 7 ^ + 7 ^ i n units of 7 ^ 0) for the 
confined modes p = 1,..., 5. Dashed lines: Landau damping. 
In the calculation d/lo = 0.25 . 

The dependence of 7H) (solid lines) for p = 1,..., 5 on 
the laser intensity is shown in Fig. [6] For sake of com¬ 
parison the Landau damping contribution is represented 
by dashed lines. In the figure it is observed that the total 
damping presents the same behavior as the Beliaev decay 
(see Fig.UJ), also, that 7^ shows resonant transitions for 
~ 140. 


IV. DISCUSSION OF THE RESULTS AND 
CONCLUSIONS 


FIG. 5. (Color online) Total damping 7 ^ = 7 ^ + in 
units of y *- 0 - 1 for the Goldstone modes p = 1, 2,..., 8 versus the 
dimensionless self-interaction parameter A. Landau damping 
for the modes are shown for equivalent color dashed lines. 

as a function of A for the first fifth excited states. The 
Landau contribution 7)) ' is represented by dashed lines. 
As in the case of the Beliaev process, yhh/yl 0 ) ^ |A| for 
small values of the self-interaction atom-atom parameter, 
while for large values of |A|, the function ylH/yl 0 ) has 
a maximum at certain A p value. We note that ^/ r P> is 
smaller than 7^ for all excited states p = 2,..., 5. For the 


From the reported calculations, two main results can 
be easily deduced: a) The behavior of the damping rates 
7 g \ 7^ and 7^ with the temperature, and b) the eval¬ 
uation of the renormalized confined phonon frequencies. 
Equations (101 and allow a good approach for all 
temperature regime. Figures [t] and ^depict 7^ and 7^ 
decays for A = 2 and Vo = 0 as a function of the reduced 
temperature ksT/hui 0, respectively. From Fig. ^ we have 

(P) 

that, in the range of temperature here considered, ^ > 
7'hh for all excited states p > 1. 


In the thermal regime where fcgT >> Hoj 0 , i.e. high 
temperature or weak quantum confinement, we obtain 
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FIG. 7. (Color online) Reduced Beliaev decay rate 7 ^/ 7 ^ 
versus ksT/hui 0 for the Goldstone modes p = 2, Open 
squares represe nt th e solution for the thermal regime as de¬ 
rived from Eq. (141. Open diamonds correspond to the low 
temperature limit according to Eq. (15 1 . 


FIG. 8 . (Color online) Reduced total damping 7 ^/ 7 ® for 
the Goldstone modes, p = 1,2,3,4 as a function of tempera¬ 
ture. Dashed lines: Landau damping. Weak, k B T » ftoJo, 
and strong confinement, k B T << hcoo, limits are represented 
by open squares and open diamonds, respectively (see text). 


that given by Eq. (|13|, is reduced to 


7 (p) = 7 (o) | A | 


k B T 

hlOn 


A « + 


s(0)> o) 4^ o> ) 


(14) 


(r) (r) 

where the Landau, Ap , and Beliaev, Bp , coefficients 


(r = 0,1 ) are temperature independent (see Eqs. ( A6) 
and (B3)). The asymptotic behaviors for 7^ and 7^ 
are displayed by open squares in Figs. [7] and [8] respec¬ 
tively. In the case of the Beliaev damping and by com¬ 
parison with the exact result in Fig. [7J we can see that 
the thermal regime supported by Eq. ( |A5[ ) , is a better 
approach for the lower excited states. For example, for 


p = 2, Eq. (A5 ) provides a good result if k B T > 0.5fiwo, 


while for p = 8, we have that the function (A5 ) is valid if 
k B T > 1.3Hojo- In contrast, for the total damping, i.e. for 


higher values of p, the limit ( 141 becomes a better ap¬ 
proach. Thus, for p = 2 and 4, Eq. (141 is good enough 
if k B T > hcoo and k B T > 0.6Hojo, respectively. Notice 


that for p = odd numbers, the thermal regime (141 is 


even better. These facts are explained by the presence of 
the Landau matrix element, AlP , in the total damping 
calculation for T 7^ 0 K. 

The linear character of the damping rate with T has 
been tested experimentally in the atomic gas of NaP In 
the experiment of Ref. |4| the condensate was loaded in a 
trap where the transversal frequency w r >> ujq. Hence, 
we can argue that we are in presence of a quasi-ID con¬ 


densate. The excitation frequency employed in the exper¬ 
iment was of uj ex = 1.58u;o and ujq = 2ttx 19.3 Hz. Follow¬ 
ing the Bogoliubov excitation spectrum of Eq. ([5]) , the ex¬ 
citation frequency w ex corresponds to the p = 2 confined 
phonon mode with a dimensional non-linear parameter 
A =3.42. Using the asymptotic expression for the high 
temperature regime, Eq. (141, we obtain 7^' = 4.4 s _1 
for T = 200 nK and 17.6 .s~ 1 for T = 800 nK which agree 
quite well with the reported experimental values of 4.4 
s -1 and 18 s . In the evaluation was assumed a con¬ 
densate of 3500 atoms and from the value of A =3.42 we 
extract an effective ID coupling constant gi = 3.7 x 10 -25 
eVm. 

At very low temperature or strong confined regime, i.e. 
k B T « hu>[) 1 from Eqs. (A7l and (B4) follows that the 


total damping of the exited mode p can be cast as 


7 (p) = 


A«(T) + i 


(15) 


Here, the coefficients and Bp 1 ' decay exponentially 
with T and yw) is almost constant independent of the 
temperature. Comparing the results of Eq. (15) with the 


theoretical calculations for 3D or 2D homogeneous sys¬ 
tems, where the trap potential and confined effect are 
neglected, we found for the Landau damping a different 
behavior. Reference [?3j reports the law 7^ ~ T 2 , while 
the limit of 7 B ~ T 4 is predicted in Ref. [ |T7, EDS and 
I20| . The quantum limit or very low temperature for Be- 
liev damping and the total decay as a function of re¬ 
duced temperature, are represented by open diamonds 
in Figs. [7] and [8j From the figures it can be notice that, 
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FIG. 9. Renormalized phonon frequencies, Re{n p }=Aoj p in 
units of 7^ as a function of the reduced self-interaction pa¬ 
rameter A for three values of the temperature ksT/hcoo = 
0,1, 2 and phonon states p = 1,2,3,4. 


In conclusion, we evaluated the damping rates of con¬ 
fined phonon modes of ID condensates in a harmonic trap 
potential loaded in an optical lattice. We remarked the 
influence of the spatial confinement potential on the col¬ 
lective oscillations and on the damping rates as a function 
of the temperature. The presence of an optical lattice as 
an external field, allows to manipulate the decay rate of 
the condensate. The damping 7^ can be turned on or 
turned off as a function on of the laser intensity. Also, 
for a given excited state p and tuning the laser intensity, 
it is possible to get a set of transitions, oj p —1- ojj ± 
reaching to a resonant effect for the total lifetime 1 /7W. 
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Appendix A: Beliaev matrix element 


for ksT « O.GHujo, the asymptotic approach given by 
Eq. (15) reproduces quite well the decay processes. 


An import result is the knowledge of the excited fre¬ 
quency shift as a function of the condensate parameters 
and the applied laser intensity. The real part of the self¬ 
energy in Eq. © allows to an analytical expression for 
the renormalized excited frequency, Re{Tr p }=Au p , as a 


function of A, Vo and T. Thus, 
Auj p = 70 |A| ^ [(1 + fi + fj) 


B, 


(u} p - Uj - U}j)ujQ | 

{oj p — u>i — Wj) 2 + e 1 

l ff f \ IA—I 2 (Wj - Wi - Wp)w 0 

2 Ui ( w _ - Lu p ) 2 + e 2 


• (16) 


Figures © shows the dependence of the dimensionless 
self-interaction parameter A on the renormalized discrete 
phonon frequencies A uj p = to — u p for the reduced tem¬ 
perature values ksT/hojo = 0, 1 and 2. We conclude that 
for the attractive regime (A < 0), the renormalized shift 
A oj p > 0, while the opposite result is obtained for the 
repulsive interaction, i.e. A u> p < 0 if A < 0. This behav¬ 
ior is understood by the dependence of A uj p in Eq. (16) 


on the label spacing A p ’^ = (oj p — 07 — uij) /wq and 

A(*’ p ^ = (1 ujj — LJi — u p ) / Uq as a function of A. According 
to the results of Appendixes A anb B, we have that in the 
thermal regime, A ui p is proportional to ksT/Hui o- Thus a 
linear increase or decrease of the excited frequency with 
the temperature is predicted for attractive or repulsive 
interaction between atoms, respectively. 


After substituting the perturbed wave function |it p ) 
and | u p ) in the matrix element (|9|, and neglecting terms 
of the second order or higher in Aand Vq , for the function 
we get 

M^(A,Vo) = + fi + ^ , (Al) 

i,3 

where B^j = Topij A Fpij H - VoBpij , Tivij ( ^ p ^ j = 
even number) is reported elsewhere 

Fpij = a pi j + fet. + bf v + 2{b pij + b~- p + b~ pi ) , (A2) 


Hpij Cpij H - dpij + dijp + djpi , (A3) 


E -^2m000-^2mpzj # 7 + 

- ; 0 , 

m^O 


2m 


= £ 


T T 
- L 00 pm - L 0 mij 




m d= p 

E dO^md^Omij m Qp^rn-^Omij 

2m ’ 2(m — p) 


m^O 


m^p 


(A4) 

with the parity condition p + i + j = even number, and 
the Lorenzian function 


£(*) = - 


U 0 £ 


(Wp =F OJi -Uj) 2 +£ 2 


In the limit of thermal regime, the probability Aip 8 ' 1 is 
reduced to 




gW + if B (°) 

p 12 \k B T) p 


(A5) 


































with 


If k B T << Hco o, the Landau probability process can be 
approached to 


4 r> = E 

i,3 


UJi + CUj 
UJo 



For low temperature we have 


0 , 1 ) . 

(A6) 


M ( p B) ( T) = + BW (T) . (A7) 


Mf\T) 


k B T 

huio 


/t (!) 


1 

12 


fwp \ 

k B Tj 


2 


j(0) 


(B2) 


where 

®p r) = X] ( ex P (-hoji/k B T) A exp (—hu}j/k B T)) r 
i,3 

x|B^| 2 4+>; (r = 0,1) . (A8) 


Appendix B: Landau matrix element 

Using the wave function | u p ), |u p ) and Eq. ( jlij ) and 
neglecting terms higher than A and Vp we have for A4p ' 

H L) = EcA-/i) ra 2 4 _) . ( B1 ) 

where A?j = Topij A Dpij A VpGpjj . 

Dpij &pij A bpij A b^p A 2(b pi j A b^j p A j p i ) , 

Gpij — Cp,y A dpij A dijp A djpi . 


with 




(* 


While, for the weak confinement, k B T » Hcop, 
tained that 


= 0 , 1 ) . 

(B3) 
it is ob- 


} (T) = A« = ]T [eM-f^i/ksT) 

ij 

— exp (—hujj/kBT)} 

x|A“| 2 4-). (B4) 
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